Stripe phases in the two-dimensional Falicov-Kimball model 
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The observation of charge stripe order in the doped nickelate and cuprate materials has motivated 
much theoretical effort to understand the underlying mechanism of the stripe phase. Numerical 
studies of the Hubbard model show two possibilities: (i) stripe order arises from a tendency toward 
phase separation and its competition with the long-range Coulomb interaction or (ii) stripe order 
inherently arises as a compromise between itinerancy and magnetic interactions. Here we determine 
the restricted phase diagrams of the two-dimensional Falicov-Kimball model and see that it displays 
rich behavior illustrating both possibilities in different regions of the phase diagram. 



71.28+d, 71.30+h, and 71.10-Hf 

The discovery of charge stripes in the nickelate pj and 
cuprate Q materials has encouraged much theoretical 
work to explain the underlying physical principles be- 
hind the stripe order and to examine whether the stripes 
are related to the mechanism for high-temperature su- 
perconductivity. There are two schools of thought for 
the physics that drives stripe formation: (i) Kivelson, 
Emery, and coworkers S propose that strongly corre- 
lated systems have a natural tendency toward phase sep- 
aration and the inhomogeneous spatial charge ordering 
arises from a competition between this tendency to phase 
separate and the long-range Coulomb interaction which 
does not allow the electron density to stray too far from 
its average; and (ii) Scalapino and White 0] propose that 
the stripe order arises from a competition between kinetic 
and exchange energies in a doped antiferromagnet which 
does not require long-range Coulomb forces to stabilize 
the stripes. Despite a large amount of numerical work 
ranging from high-temperature expansions jq| , to Monte 
Carlo simulations || , to exact diagonalization B , as well 
as semiclassical Hartree-Fock theory M , no consensus has 
been reached about the region of stability for the phase 
separated states or the mechanism for stripe formation. 

Here we take an alternate point of view. Rather than 
try to prove phase separation in the Hubbard or t — J 
models, we choose an even simpler model — the spinless 
Falicov-Kimball model |)| , which can be analyzed exactly. 
The relation of the Falicov-Kimball model to the Hub- 
bard model is analogous to the relation between the Ising 
and the Heisenberg models of magnetism (the Falicov- 
Kimball model can be viewed as a Hubbard model where 
the down-spin electrons are frozen and do not hop) . The 
Hamiltonian is 



H 



etc, 



U V \c\c x w x , 



(1) 



with c x (c x ) the creation (annihilation) operator for a 
spinless electron at site x and w x is a classical variable 
that denotes the presence (absence) of an ion at site x 



when it is equal to 1 (0), respectively. The hopping oc- 
curs between nearest-neigbors on a square lattice and the 
interaction strength is denoted by U . For any given con- 
figuration of ions {w x } the ground state for N e electrons 
is determined by diagonalizing a one-body operator given 
by the above Hamiltonian, and filling in the lowest N e 
states. We typically are interested in the ground-state 
configuration of the ions for a given number of electrons 



4, 



and a given number of ions Ni = J2 X 



(or their densities p e — N e /N and pi — Ni/N respec- 
tively) . 

The possibility of phase separation when U — > oo was 
proposed in 1990 |10| and is called the segregation prin- 
ciple. It was proved in the one-dimensional case flTT| ], 
in the infinite-dimensional case |l2| , and recently in the 
general case Jj"3| . This phase separation is a special type 
of phase separation, often referred to as the segregated 
phase, where the electrons and ions avoid each other, and 
reside in separate domains. It is the analog of the ferro- 
magnetic state in the Hubbard model. At half filling for 
the electrons and the ions, the ground state is known to 
be the chessboard phase |l4],[t5| ■ This state is the analog 
of the antiferromagnetic state in the Hubbard model. 

The question we pose is, if we fix the number of elec- 
trons to equal the number of ions, then what are the 
stable phases as a function of the number of electrons. 
This problem is the analog of doping the S z — phase 
of the Hubbard model away from half filling. Since we 
know the system phase separates as p e — > and is in the 
chessboard phase for p e — > 1/2, we have the interesting 
situation of determining how the transition is made from 
the segregated phase to the chessboard phase, which are 
two phases that are about as different from each other as 
possible. 

It is well known that the ground-state phase diagram 
of the Falicov-Kimball model typically includes a large 
number of different phases as functions of the particle 
concentrations and the interaction strength. As such, it 
seems unlikely that one can rigorously analyze ground- 
state phases in the general case, although a number of 
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interesting results exist for special cases |T(| (usually for 
large U) . Instead, we work with the restricted phase di- 
agram technique, where we consider all possible periodic 
phases for which the number of sites per unit cell -/Vo is 
less than or equal to N c = 16. The technique first identi- 
fies all nonequivalent periodic phases, which number 23 
755 in our case. Then for each periodic phase in our trial 
set, we calculate the total energy at a given value of U. 



Since, unlike the ID case (see 17|), there are no exact for- 
mulae known for the density of states of general periodic 
phases in 2D, we have performed a numerical solution of 
the corresponding eigenvalue problem to determine the 
bandstructure (see Ref. [|l8| for the details) . This involves 
finding the eigenvalues of an No — dimensional matrix 
for each value of k in a two-dimensional grid covering the 
Brillouin zone. Our calculations were performed with a 
k — space grid of 110x110 points for each phase. Then 
the grand canonical phase diagram is constructed as a 
function of the electron and ion chemical potentials, and 
finally the grand-canonical phase diagram is translated 
to the canonical phase diagram for arbitrary pi and p e . 
This procedure assures thermodynamical stability of all 
phases (both periodic and their mixtures) present in the 
resulting canonical phase diagram. The stability problem 
is discussed extensively in Ref. jl9| , where the canonical 
phase diagrams of the ID Falicov-Kimball model were 
studied. Finally, we make the restriction pi = p e in the 
canonical phase diagram to analyze the problem at hand. 

The ground-state phase diagram is quite complex. We 
find many different stable phases occur, which can be 
classified into a number of different categories: (i) the 
empty lattice (pi — and p e ^ 0) denoted E; (ii) the full 
lattice (pi = 1 and p e = 0) denoted F; (iii) the chessboard 
phase (pi = p e = 1/2 and ions occupy the A sublattice 
only) denoted Ch; (iv) diagonal neutral (fli = 1 — p e ) 
stripe phases (the ions are arranged as diagonal chess- 
board phases separated by fully occupied striped regions 
with a slope of one, or equivalently, Ch phases separated 
by diagonal antiphase boundaries) denoted DNS; (v) di- 
agonal non-neutral {pi ^ 1 — p e ) stripe phases (the ions 
are arranged as diagonal chessboard phases but separated 
by empty striped regions with a slope of one) denoted 
DS; (vi) axial non-neutral stripes (pi 1 — p e and ions 
arranged in stripes parallel to the y-axis and translation- 
ally invariant along the axis) denoted AS; (vii) axial non- 
neutral chessboard stripes (pi =/= I — p e and ions arranged 
in stripes of the chessboard phase oriented parallel to 
the y-axis) denoted AChS; (viii) other neutral phases 
(pi = 1 — p e but the arrangement is not in any sim- 
ple stripe-like phase) denoted N; (ix) non-neutral four- 
molecule phases (where pi =/= I — p e and empty sites are 
arranged out of "bound" four-molecule squares) denoted 
4M; and (x) truly two-dimensional non-neutral arrange- 
ments of ions (where pi ^ l—p e and the ions are arranged 
in a fashion that is not stripe-like, but rather requires 
a two-dimensional unit cell to describe it) denoted 2D. 



Generically, the phase diagram includes mixtures of two 
or three of the periodic phases (iii-x), or of one or two 
periodic phases and the empty lattice (i). The empty 
lattice is usually needed in the phase mixtures to ensure 
that the average electron and ion fillings are equal for 
the mixtures. For occasional values of the filling, we do 
find single periodic phases to be stable, but this feature is 
rare. When the filling is sufficiently far from half filling, 
the system is in the segregated phase, which is a mixture 
of the empty lattice (pi = 0, p e ^ 0) and of the full lattice 

(fli = 1, p e = 0). 

The number of phases stable or appearing in mixtures 
tends to grow as U is decreased in magnitude. For U = 8 
there are 25 phases: E (1); F(l); Ch (1); DNS (6); and 
N (16). For U = 6 there are 30 phases: E (1); F (1); Ch 
(1); AS (20); N (6); and 2D (1). For U = 4 there are 42 
phases: E (1); F (1); Ch (1); AS (35); N (1); and 2D (3). 
For U = 2 there are 38 phases: E (1); F (1); Ch (1); AS 
(14); DNS (4); 4M (2); and 2D (15). For U = 1 there are 
50 phases: E (1); F (1); Ch (1); DS (1); AS (9); AChS 
(6); 4M (3); and 2D (28). 

(a): 2 x 1 (1/2) (b): 9 x 1 (5/9) (c): 15x1 (7/15) 



(d):3x1 (2/3) (e):7x2 



(3/7) (f):9x1 (7/9) 



(g): 15x1 (11/15) (h): 3x3 (4/9) (i): 3x2 (1/6) 
• • • ■ • • • 



FIG. 1. Characteristic configurations representing phases 
displayed in the phase diagram: (a) Ch, (b) DNS, (c) DS, (d) 
AS, (e) AChS, (f) N, (g) 4M, (h) 2D, and (%) another 2D. The 
large dots correspond to ion-occupied sites and the small dots 
correspond to ion-vacant sites. The numbers written above 
each configuration indicate the dimensions of the unit cell; 
the corresponding ion densities are given in parentheses. 

The empty lattice, the full lattice and the chessboard 
phase are present in the phase diagram for all U; the 
diagonal neutral stripes generally for large U (U > 7) 
(occasionally they can appear for moderate U as well — 
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as occurs for U = 2); the axial stripes for moderate and 
small U (U < 7); the axial chessboard stripes for small U 
(U < 3); the other neutral phases for large U (U > 4); the 
four-molecule phases for small U (U < 3); and the two- 
dimensional phases for moderate and small U (U < 7) 
(growing significantly in number as U decreases). The 
total number of phases appearing in the phase diagram 
are too large to illustrate here. Instead, we choose a 
number of illustrative examples of each of the different 
types. 

The empty lattice (E) corresponds to a lattice with 
no ions and the full lattice (F) corresponds to a lat- 
tice fully occupied by ions. The chessboard phase (Ch) 
corresponds to the case where ions (large dots) occupy 
the A sublattice and empty sites (small dots) occupy the 
B sublattice as shown in Fig. 1(a). A diagonal neutral 
stripe phase (DNS) is shown in Fig. 1(b). It is typical 
of the phases we see with regions of chessboard phase 
separated by antiphase boundaries. This configuration 
has pi = 5/9 and p e — 4/9 and every ninth diagonal 
is an antiphase boundary. An example of a nonneutral 
diagonal stripe phase (DS) which has pi = 7/15 and 
p e = 7/15 is shown in Fig. 1(c). Like the neutral stripes, 
these have antiferromagnetic regions but they are sepa- 
rated by diagonal empty (and not full ) lattice stripes; 
the empty lattice stripes are width two now. An ax- 
ial stripe phase (AS) is shown in Fig. 1(d). It is typ- 
ical of these phases with fully occupied (ferromagentic) 
stripes separated by the empty lattice. This configura- 
tion has pi = 2/3 and 0.189 < p e < 0.283; other sta- 
ble phases have wider bands of occupied stripes, which 
are always separated by an empty lattice stripe of width 
one. An example of an axial chessboard stripe phase 
(AChS) is shown in Fig. 1(e). These phases consist 
of regions of the chessboard phase separated by width 
one empty lattice stripes. The width of the chessboard 
phase ranges all the way down to two. This particular 
phase has pi — 3/7 and p e = 3/7. An example of an 
other neutral phase (N) is depicted in Fig. 1(f). These 
phases can be viewed as striped phases, alternating oc- 
cupied (ferromagnetic) and empty lattice stripes, with 
the stripes having slope different from 0, 1, or oo. In 
this example, pi = 7/9 and p e = 2/9; the slope of the 
empty lattice stripe is 1/2. A four-molecule phase (4M) 
is shown in Fig. 1(g). These phases are two-dimensional 
tilings of four-molecule empty-site squares inside an oc- 
cupied "latticework". This phase has pi = 11/15 and 
p e = 1/15. The truly two-dimensional phases (2D) are 
more difficult to classify — some, like in Fig. 1(h), have 
no stripe-like character at all (pi =4/9 and p e = 4/9), 
others can be described with a stripe-like picture such as 
in Fig. l(i) (p t = 1/6 and 0.32 < p e < 0.34) which can 
be viewed either as slope 2/3 ferromagnetic stripes in an 
empty lattice, or as antiferromagnetic axial stripes, etc.; 
in the end we classify these phases as two-dimensional. 
The summary of all of the stable phases will appear in a 



longer publication j2C§ 
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FIG. 2. Phase diagram of the 2D FKM along the p e = pi 
line for a set of U values. The dot-dashed line separates the 
segregated phase and mixtures of the E phase with periodic 
phases. Below this line only the E and F phases coexist (in 
fact the F phase is also stable slightly above the line, forming 
three-component mixtures within narrow intervals of densi- 
ties). Vertical line segments mark intervals of the densities 
where the corresponding phases (of a given class) are involved 
in the formation of stable phases. (All the segments, except for 
those corresponding to the E phase, are moved slightly right 
or left of the value of U/(U + 8).) Bold lines correspond to 
various types of stripe phases. If only one category is marked 
within an interval, only phases belonging to this category can 
form stable phases ( either single periodic phases or mixtures 
of two or three phases from within the same class). If seg- 
ments corresponding to different categories overlap within a 
given interval, then only mixtures of phases belonging to those 
categories can form stable phases. 

The phase diagram is shown in Fig. 2 and its simpli- 
fied, schematic version, in Fig. 3. The phase boundary 
between the segregated phase and mixtures with peri- 
odic phases, approaches half filling as U gets large, as 
expected. When U = 8, we find the mixtures, when 
doped just away from half filling, occur between the 
chessboard phase, the empty lattice and diagonal neu- 
tral stripe phases. This picture is similar to those that 
propose the phase separation scenario, but there is no re- 
quirement of the long-range Coulomb interaction to gen- 
erate the stripes — they also occur as part of the periodic 
phases that compose the different stable mixtures. As 
the system is doped further away from half filling, mix- 
tures with other neutral phases occur, before the system 
fully phase separates. 

As U is reduced, the phase diagram becomes more 
complicated. Near half filling, the chessboard phase is 
always one of the phases in the stable mixtures, but we 
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find the empty lattice and diagonal neutral phases dis- 
appear and are replaced by other neutral or 2D phases 
in the mixtures. Then as U is reduced further, the diag- 
onal neutral stripe phases occasionally re-enter into the 
mix, replacing the N or 2D phases, as can be seen for 
U = 2. Finally, for smaller U, the mixtures are be- 
tween the diagonal non-neutral stripes and the chess- 
board phase only. Farther away from half filling, the 
behavior is even more complex, with axial non-neutral 
stripes first entering near the segregation boundary, and 
then the 4M phases appearing as U is reduced further. 
The axial neutral chessboard phases also appear for small 
U. It is in this moderate-to-small U region where we see 
many stripe phases form due to a delicate balance be- 
tween kinetic-energy and potential-energy effects. This 
is the alternate scenario for stripe formation, that does 
not require phase separation or the long-range Coulomb 
interaction. 
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FIG. 3. Schematic phase diagram with highlighted regions 
of stability of the stripe phases. 

In summary, we have calculated the restricted phase 
diagram of the two-dimensional Falicov-Kimball model. 
We constrained the system to have the same number of 
localized and itinerant particles, which is the analog of 
the S z = state of the Hubbard model. We find that, in 
addition to the phase separation of the segregated phase, 
the system generically forms a number of different stripe 
phases. For large U, we find the stripes to be diagonal 
stripes and appearing only close to half filling. As U 
is reduced, we find diagonal stripes changing into axial 
chessboard stripes and then axial stripes as the system is 
doped further away from half filling. In addition, we find 
a number of truly two-dimensional phases present as well. 
While we cannot say anything about what happens in the 
Hubbard model itself, our results suggest that one should 
expect a complex phase diagram when stripe phases are 
present and see a competition between the stability of 
the stripes and other, more two-dimensional, structures. 



Of course, we can only speculate on the behavior of the 
Hubbard model, but this exact solution of the Falicov- 
Kimball model shows that stripe formation is generically 
a complicated occurrence, that can have a number of 
competing states close in energy. We do expect, however, 
that the phase diagram simplifies for the Hubbard model, 
since the system can quantum-mechanically fluctuate be- 
tween different low-energy "ionic" configurations. 

Acknowledgements: R.L. and G.B. acknowledge sup- 
port from the Polish State Committee for Scientific Re- 
search (KBN) under Grant No. 2P03B 131 19 and 
J.K.F. acknowledges support from the NSF under grant 
No. DMR-9973225. We also acknowledge support from 
Georgetown University for a travel grant in the fall of 
2001. 



[1] C.H. Chen, S.-W. Cheong, and A.S. Cooper, Phys. Rev. 

Lett. -71, 2461 (1993); J.M. Tranquada, et al., ibid. 73, 

1003 (1994); Phys. Rev. B 52, 3581 (1995); V. Sachan, 

et al., ibid. 51, 12742 (1995). 
[2] J.M. Tranquada, et al., Nature (London), 375, 561 

(1995); Phys. Rev. B 54, 7489 (1996); Phys. Rev. Lett. 

78, 338 (1997); H.A. Mook, P. Dai, and F. Dogan, Phys. 

Rev. Lett. 88, 097004 (2002). 
[3] V.J. Emery, S.A. Kivelson, and H.Q. Lin, Phys. Rev. 

Lett. 64, 475 (1990); E.W. Carlson, et al., Phys. Rev. B 

57, 14704 (1998); L.P Pryadko, S.A. Kivelson, and D.W. 

Hone, Phys. Rev. Lett. 80, 5651 (1998); L.P. Pryadko, 

et al., Phys. Rev. B 60, 7541 (1999). 
[4] S.R. White and D.J. Scalapino, Phys. Rev. Lett. 80, 1272 

(1998) ; ibid., 81, 3227 (1998); Phys. Rev. B 60, R753 

(1999) ; ibid., 61, 6320 (2000). 

[5] W.O. Putikka, M.U. Luchini, and T.M. Rice, Phys. Rev. 

Lett. 68, 538 (1992); W.O. Putikka and M.U. Luchini, 

Phys. Rev. B 62, 1684 (2000). 
[6] A.C. Cosentini, et al. Phys. Rev. B 58, R14685 (1998); 

M. Calandra, F. Becca, and S. Sorella, Phys. Rev. Lett. 

81, 5185 (1998); C.S. Hellberg and E. Manousakis, Phys. 

Rev. B 61, 11787 (2000). 
[7] E. Dagotto, Rev. Mod. Phys. 66, 763 (1994); C.S. Hell- 
berg and E. Manousakis, Phys. Rev. Lett. 83, 132 (1999). 
[8] J. Zaanen and A.M. Oles, Ann. Physik 5, 224 (1996); D. 

Cora, K. Rosciszewski and A.M. Oles, Phys. Rev. B 60, 

7429 (1999). 

[9] L.M. Falicov and J.C. Kimball, Phys. Rev. Lett. 22, 997 
(1969). 

[10] J.K. Frecricks and L.M. Falicov, Phys. Rev. B 41, 2163 
(1990). 

[11] P. Lemberger, J. Phys. A 25, 715 (1992). 

[12] J.K. Frecricks, Ch. Gruber, and N. Maoris, Phys. Rev. B 
60, 1617 (1999). 

[13] J.K. Frecricks, E.H. Lieb, and D. Ueltschi, Phys. Rev. 
Lett. 88, 106401 (2002); Commun. Math. Phys. (to ap- 
pear) (2002). 



4 



[14] T. Kennedy and E.H. Lieb, Physica 138A, 320 (1986); 

ibid 140A, 240 (1986). 
[15] U. Brandt and R. Schmidt, Z. Phys. B 63, 45 (1986); 

ibid. 67, 43 (1987). 
[16] T. Kennedy, Rev. Math. Phys. 6, 901 (1994); T. 

Kennedy, J. Stat. Phys. 91, 829 (1998); K. Haller, PhD 

dissertation, Univ. of Arizona (1998); K. Haller, Com- 

mun. Math. Phys 210, 703 (2000); K. Haller and T. 

Kennedy, J. Stat. Phys. 102, 15 (2001). 
[17] R. Lyzwa, Phys. Lett. A 164, 323 (1992); Physica A 192, 

231 (1993). 

[18] G. I. Watson and R. Lemahski, J. Phys.: Condens. Mat- 
ter 7, 9521 (1995). 

[19] Z. Gajek, J. J§drzejewski and R. Lemahski, Physica A 
223, 175 (1996). 

[20] R. Lemahski and J.K. Freericks, unpublished. 



5 



